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COUPLED ROTOR AND FUSELAGE EQUATIONS OF MOTION 


William Warmbrodt 
Ames Research Center 

SUMMARY 

The governing equations of motion of a helicopter rotor coupled to a rigid 
body fuselage are derived. A consistent formulation is used to derive non- 
linear periodic coefficient equations of motion which can be used to study 
coupled rotor/fuselage dynamics in forward flight. Rotor/f: selage coupling 
is documented and the importance of an ordering scheme in devrivlng nonlinear 
equations of motion is reviewed. The nature of the final equations and the 
use of multiblade coordinates are briefly discussed. 

I. INTRODUCTION 

Helicopter rotors are complex aeroelastic systems where structural, 
inertia and aerodynamic forces interact to determine rotor behavior. To 
mathematically model this problem, a consistent formulation of the governing 
equations of motion must be used. Additionally, recent research has indicated 
that rotary wing aeroelasticity is inherently nonlinear} Consequently this 
aeroelastic problem requires a careful derivation of the dynamic equations 
of equilibrium so that moderate blade deflections based upon the assumption 
of small strains and finite slopes are properly incorporated in the rotor 
math model. 

Mathematical models of single isolated helicopter blades are well 
2 3 

developed.’ The results of these investigations have shown that the effect 
of nonlinear terms on blade aeroelastic stability is important. However, 


there has been relatively little emphasis on developing consistenUy derived 
nonlinear math models detailing coupled rotor/fuselage dynamics. 

Considerable work has been done to analytically model coupled rotor/ 

fuselage systems. These models include large general purpose helicopter 

, 4-7 

analyses and math models which study aeromechanical phenomenon for a 

8—10 

particular rotor configuration . The first group of analyses utilize 
a general rotor model and system degrees of freedom. Reference 11 presents 
a review of such comprehensive helicopter analyses and the level of technology 
used in the math modeling procedure. Despite the sophistication of these 
models, there has been little effort to derive consistent nonlinear equa- 
tions of motion. The second group of analyses limit attention to a given 
rotor configuration and flight regime or restrict fuselage motion. 

This report will present the derivation of a set of equations which 
model the basic dynamics of a coupled hingeless rotor/fuselage system. The 
math model is not intended to represent all the very complicated interactions 
the aircraft experiences in flight. Rather, this study is an attempt to 
mathematically model the complex dynamic interaction between the rotor and 
fuselage. This is performed by deriving a set of nonlinear equations of 
motion while clearly documenting the rotor/fuselage matching procedure. 


II. FORMULATION OP THE PROBLEM 

A. General 

The configuration considered in this study is shown in Figs. 1 and 2. 
The rotor has N blades with distributed mass and stiffness properties. 
Each blade can have out-of-plane bending, , inplane bending, ^ and 


2 


torsion, » deflections. Rotor geometry Is shovm in Figs, 3 and 4. 

The helicopter fuselage is modeled as a rigid body with three transla- 
tional and three rotational degrees of freedom. The fuselage center of mass 
is assumed to be located in the vertical centerline plane of the fuselage. 

B. Basic Assumptions 

The following assumptions are used in deriving the equations of motion. 

1) The blades are connected to the hub at an offset distance e from 
the axis of rotation. 

2) The blade feathering axis is preconed by an angle The blade 
has no torque offset and no angle of sweep or droop. 

3) The blade cross section is symmetric about the major principal 
axis. 

4) The blade cross section can have four distinct points: center 

of gravity, aerodynamic center, tension center, and elastic 
center. The elastic axis is assumed to be a straight line and 
coincident with the feathering axis, 

5) The blades can have distributed pretwist about the undeformed 
elastic axis. 

6) The blades are assumed to be composed of an isotropic linear 
elastic homogeneous material. 

7) Each blade can bend in two mutually perpendicular directions 
normal to the elastic axis and torsionally twist. Moderate 
deflections are assumed resulting in small strains and finite 
rotations. 


3 


8) During deformation, the blade cross sections reronin planci and 
normal to the elastic axis (Euler-BerrAOulli hypothesis). 

9) Two-dimensional quasi-steady aerodynamic loads are used. Apparent 
mass, compressibility, and stall effects art ittcluded; dynamic 
stall effects are neglected. 

10) The rotor is assumed to operate at a constant sp.ied of rotation. 

11) All structural damping forces are assumed to be of a viscous type. 

12) The rotor shaft is assumed rigid. No control system tlexibility 
or engine dynamics are modeled. 

13) The fuselage is modeled as a rigid body with total mass M and 

F 

momf'Uts of Inertia I„, I , I and product of inertia I 

A A Z Za 

“ ^XY “ * 

14) The fuselage center of mass is located in the vertical centerline 
plane of the fuselage. 

15) The fuselage has translational degrees of freedom R^, R^, in 
three mutually perpendicular ddrections and »*hree rotational degrees 
of freedom 0„, 0„, 0„ . 

16) The helicopter is assumed to be in equilibrium forward flight. 
Additional assumptions used in the derivation are given in the text. 

C. Ordering Scheme 

lu the process of deriving the set of governing nonlinear partial dif- 
ferntial equations, a considerable number of terms which may be small and 
negligible are encountered. To neglect higher order terms in a systematic 
manner, an ordering scheme is used^. All quantities are assigned an order 
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of magnitude in corws of a nondimeneional quantity, t, which represents 

typical elastic blade slopes. The ordering scheme is then used with the 

2 

assumption tliat terras of 0 (e ) are negligible when compared to terras of 

2 

order one, i.e, 1 + 0(e ) ■ 1, 

When raeieling the dynamics of coupled rotor/support systems, the 

magnitude assignment of fuselage motion is very important. The assigned 

order of magnitude of fuselage motions must allow the moth model to adequately 

represent actual system behavior. However, the ordering scheme should bo 

such to restrict the total number of terms appearing in the final equations 

to those which are important. This will not only allow the solution of the 

final equations to be less cumbersome, but also yield greater physical 

insight into the nature of the problem. Consequently, it is assumed in this 

study that fuselage rotations and perturbation translations are of order 
3/2 

e . This will allow the model to be valid for a helicopter in equilibrium 

flight undergoing small fuselage motions. Moderate blade deformation is 

retained. Additionally, the analysis will remain general to adequately 

model the complex dynamics between the rotor and the fuselage. On the other 

hand, this assignment restricts Che tremendous number of fuselage-dependent 

rotor forcing terras. If fuselage rotations and nondimensional displacements 

are assigned a lower order of magnitude, i.e., 0(e), the expressions for 

the rotor/fuselage matching conditions are practically unmanageable. The 

3/2 

present assignment of 0(e ) will restrict the final equations of motion 

Co a length where physical understanding is still possible. 
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The uaalgned orders of magnitude of the quantities used in this study 


are 
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The systaraatxu applioation of this ordering scheme in the derivation procedure 
will yield a consistent set of nonlinear equations of motion. 


Ill, DERIVATION OF GOVERNING EQUATIONS 


A, General 

The math model is composed of an N-bladed elastic rotor and a rigid 
body fuselage. Each blade will have governing partial differential equations 
of motion for the flap, lag, and torsion deflections. The fuselage motion 
will be governed by the three force and three moment equations of system 
equilibrium. 

The math model of this study considers each rotor blade individually 
and the equations for the flap, lead-lag, and torsion of each blade are 
derived. Rigid body fuselage motions are included in Che calculation of 
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the bXade loadst Consequently, the equations of motion for each rotor blade 
are coupled to the fuselage through the forcing terras in the loading expres- 
sions. 

Overall system equilibrium is enforced by requiring force and moment 
balance at the rotor hub. These six equations of equilibrium are the 
governing fuselage equations of motion. 

B. Coordinate Systems 

The coordinate systems used in the derivation procedure are shown in 

Figs. 1-7. A brief description of these coordinate systems is given in 

Table 1. The Euler transformation matrix between each system is given in 

Appendix A. The unit coordinate vectors of Che jth system in the axe 

J J J 

directions are given by e , e. , e , respectively. 

jx j I jx 

C. Blade Equations of Motion 

The governing equations of motion for the kth blade of an N-bladed 
rotor r.'.i*e taken from Ref. 12. It is important to note that Ref. 12 neglects 
the influence of axial forces on the torsional rigidity of the rotor blade 
as being very small (compared to the bending stiffness). In addition, for 
rotor blades where cross sectional warping due to torsion may be significant, 
a more appropriate elasticity theory should be used to derive the blade math 
model. For the present analysis, this model is assumed to be sufficiently 
accurate. 
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The kth blade equations of equilibrium written in the undeformed coordinate 
directions, under the assumptions of Sec. II. B, are given by 

£ ('Z 22 C 05 Jyy. 5-/>7 

• sjn^Gi, cos- 6^/^) - xt-Ccc9s(%^- 

% 3Xkh n ~ ^s'zk^x ~p3Yk "" ^ 

“ t ^ Z - <^Gk ^Gk (^k^xx ^<Pk^k 'X X*) 

* B.ix^^Sin'^d-Gu ^ lyyCos'^O-G,^) 

“ Tk^r ^Gk j 

^ ^ ki ' y ^ SXk^t X. ^ 3Y k • K '7^3'Z. k. ^ ^ - 3 ^ 


4- x-7-( W^, ^^003^0 1 sind-QO 

•*• <^3Xk ^k'X^3Ek== ^ 

( 4 ) 


These represent the equations of equilibrium for the kth blade axial direction 
and the lead-lag, flap and torsion deflections, respectively. 

The distributed loading vectors which appear in Eqs. (l)-(4) may be 
written as the sum of inertia, aerodynamic and damping contributions. 


~p3k " 

Tsik 

* -Pa/Ki. * 

~PBX>k 

(5-) 


✓V' 


/V' 


%3k. ^ 

^3Xk 

r>y 

<^3Ak ^ 


w 

The detailed derivations of 

!sk 

presented below. 



C.l. Inertia Loads 

The distributed inertia loads are determined by calculating the 
inertial acceleration of a general point P on the deformed blade's cross 
section 


where the position vector r„. is written with respect to the un- 

tc 

disturbed hub location. 

The position vector is obtained by first considering fuselage 
motion. The fuselage center of mass undergoes small perturbation 
translational motion given by 

Additionally, the fuselage is assumed to undergo rigid body ro- 
tation by experiencing a yaw-pitch-roll sequence of Euler angle 
rotations, 6^, respectively. 

Fuselage motion produces hub displacement r, . As shown 
in Fig. 1, the hub is assumed to be located from the fuselage 
center of mass by 


r-H 

A 

~ Xp e. 

FX "^F^FZ 

Consequently as the fuselage experiences rigid body motion, the 
hub is displaced by 


Hp C Oy + ©X ^ RY 


1 “ 

^ F ^ ~ ^Y ^ RY 

i 

[ 

^ F ^ Y } ^ R-S 

=; 


'^Y^RV 
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Hence the position vector r , , Eq. (7), written from the un- 
disturbed hub position to the point P on the deformed blade cross 
section, is 


'"PU “ ’”V^RY ^ 

'^R‘^3Zk y-OU ^4-Yk 


^Ok ^ 4 Zk 


Oi) 


Rewriting Eq. (11) in the 2k-coordinate system using the transfor- 
mations of Appendix A yields 

’T 2 Pk =[ + r-ySin^^^ + e 4- x^, 

+ 6C^-/SpWk 

^Ok C“'^k.«~/3p -t- 2Xk 

4- { r^cos^(^ - ^ou 

^ Ok ^ ^ k . X ^ 2 Yk 

^ ^k ^ /^P (Xk ^k'^ 

^Ok (“l^p'^k'x” '^k ’ X ^k P '^k'^ 

+ ^Ok ^ ^ “ ^k/x/®P + ‘?^k^k-y /^P^}^22k 
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The angular velocity vector u , Eq, (7), Includes the fuselage 
rigid body rotations plus rotor rotation. 






d: 




A 

•h JT2- 6 . 


C13) 


Rewriting Eq. (13) in the 2k-coordlnate system yields 

*■ &y -t- e2xk 

<^X^Z ^2Vk 

+ [ 4- e--2> a4) 


Substituting Eqs. (12) and (14) into Eq. (7), performing the 
differentiations and the vector cross products, the inertial ac- 
celeration of a general point P on the deformed blade cross section 
expressed in the 2k-coordinate system is obtained, Rewriting 

A/ 

in the 3k-coordinate system yields a . , the inertial acceleration 
vector in the undeformed blade coordinate directions. Applying 
D'Alembert's principle and integrating the differential force 
acting on a differential element of mass over the blades cross 
section yields P„^, • 

.j XK 
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-P 3 Tk^-^S^.'^^ 3 Pk°>^ ' CJS) 

Introducing the identities of Appendix B, Eqs. and 

applying the ordering scheme, the final distributed inertia force 
vector of the kth blade is given by 

Tartc " - rySinl'i, * 

fCe + Xt^') -h 2 (^2 ^ ^ } ^3X \< 

■»■ cc?s+i^[-rY + 2w^&x'*' 

+ Sin4")^[ Tx 2 + dy C 4^pC e.+XUi>)] 

• * * 

- Vj^ - 2 Ui^ 4- 2 j 3 p\^^ - Ce.+ +■ v^. 

+ 2 ^ ^3Yk 

■#• C05 “ 2 V>J^ - V|^ Oy^ ~~ 2 /^P^x 

~ (:2ay^eYKe.-»-><^^] -b sin4^|^[-2eYV^, 

-v^^y 4- + /SpPy " C 2 ^y SyKe+Xj^')] 

-r^ ” ~ 2v^^^p - 2/5p ^2 Ce + xi^"^ 

-/SpCe.+.x^') - xi u}^3Ek 

" “p3lKk^3XI< 'PaxYk ^3YI< 7^3XZl<'2-32k 

06 ) 
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The distributed inertia moment is determined by taking the 
moment of the differential forces acting on *he blade cross section 
about the blade cross section elastic center and integrating over 


the blade cross section. 

Jv ^ 

Evaluating Eq. (17), introducing the idencicies of Appendix B, 
Eqs, (B-8)~(B“15) , and applying the ordering scheme yields Che 
final distributed inertia moment vector of the kch blade. 
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^VxCcse-s^j: 
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G.2 Aerodynamic Loads 


The determination of the aerodynamic loads on the rotor 
requires the resultant relative velocity between the blade and 
air. This may be written as 


V, 


V 


AU 




Ci9) 


where V., is the free stream air flow and v . is the velocity of 
the deformed blade's cross section elastic center. The free stream 
air flow includes the effect of forward flight and the induced inflow, 
which has radial and azimuthal dependence. The vector is 
given by 


C2 0) 

Rewriting Eq. (20) in the deformed blade's 5k-coordinate directions 
yields V^3j^. 

The velocity of the deformed blade's cross section elastic 
center, is given by 

a/ 


ECU 


- T£CU ^ ^ 


*^ECk'^ 


( 21 ) 


where Is given by Eq, (12) with • 0 and Wj^ is 

given by Eq# (14). Evaluating Eq. (20), rewriting in the 5k- 
coordlnate system, and substituting it along with into Eq. (19) 
yields the resultant relative velocity V . Since two-dimensional 
aerodynamic theory is used, only the velocity components in the 
plane of the deformed blade cross section are required, The final 
expressions arc 


- ^ ’ho.n of'fz ^x\ E ^ 

* ^h] “ ■" ^2 J ^22) 

+ Ca + xO&Y'*' 

+ ?y +/Sp‘^] +/^^v-A|^-r2 

- W k. ^ Yk Cy/jj, , ^ 4 - X k'^ V'k^x'^kM^ 

(23) 
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The airfoil configuration for the calculation of the steody 
aerodynamic loads is shown in Fig. 8. The steady aerodynamic lift, 


drag, and 

moment 

, per unit 

length acting on 

the kth blade are given 

by 





tsu " 




(24) 



b Cj (c^ 


(25) 


Xa 

Lsk 






(2fe) 


where 


(^Rk = 

VA ‘ ® 

(27) 


^Gk i’k 

^ysYk) 


Jjv 



(28) 

= 



(2<=f) 


The unsteady lift and moment, L^g and , are obtained from 

two-dimensional incompressible thin airfoil theory, Ref. 13. The 

assumption of quasi-steady airloads, G(k) » 1, is made. The 

influence of unsteady aerodynamics on totor blade behavior was 

first studied in Refs. 14 and 15. It was shown that a quasi-steady 

airload assumption is usually conservative when determining 

rotor blade stability. The airfoil environment for the unsteady 

aerodynamic load calculation is shown in Fig. 9. The unsteady 

lift. L , and moment, 
xxri., USk 
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, are given by 


^USk ® 1 ~ *2 ^ ^Pk 

^ ^^/ek ‘^Pk ^ ^A*^] 

^USk “ { 5 ^^A"* '1'^ *^Nk t ^Pk 

-XAb4. ^^)+ U^u «^pk[xAb- XA^^-C-l-f]} 

(3i) 


Since the unsteady lift and moment are determined using a linear 
theory, only linear terms in perturbation motion are retained in 
the unsteady load expression. Consequently, L^gj^ ^lygj^ are 
evaluated using 


^hJk^ys-Zk *" ^Vk ^ ^Sk ’*■ 

= ^cc>s- £-/4 (/Sp+ -A C&+’X^,')^yl 

(3 a) 

«Rt = - Vs-vk '! ^ 

■= ,X2l'E. jLL^fr} 4" ( ■hXi^')!^ (33) 

O'PI^ = 6^^005 4^)^ 4- 0yS^nH'|^ 4- 4- 4>J^ 

*• Pp") C3+) 
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Note that although a linear (valid up to 0 (e)) unsteady aero- 

3/2 

dynamic theory is used, terras linear in body motion ( 0(e )) 

arc retained for completeness. 


The resulting distributed aerodynamic loads written in the 
5k-coordinate system are given by 


TAS-k - f “ 

- Ds^sinC 

= [-LskC^O - 

^SYk 

^ { L +* L-uskl^szk 

‘3"A5k " 1 ^ ^USkl 

A 

~ ^AS-XU ^S~Xk 


(35) 


(3<o) 


Rewriting Eqs. (35) and (36) in the undeformed blade coordinate 
directions and neglecting the axial force component yields 




VSYU 

" -Dsk} ^3Yk 

\ ^ Sk ^U'S k } ^ 3Z k 


(37) 


^A^X.k ^A^Xk '^k'x. 3Vk 

^ 3 -Hk 


Substituting the velocity expressions into the steady and 
unsteady load expressions, evaluating, p^^j^ and 

/w V 

applying the ordering scheme yields 
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X ^ <. C-OS yOL ^ /P^ C fS p+ 
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Several comments should be made concerning the aerodynamic 

loads presented in Eqs, (39)-(42). Apparent mass terras, stall, 

and compressibility effects are included in the analysis. The 

aeroelastic phenomenon of dynamic stall, although neglected in 

the present analysis, may be incorporated in a semi-empirical 

manner. Including these effects will yield a more I'ealistic rotor 

model. Several such approaches are available in the literature 
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for aeroelastic response calculations * . Also the aerodynamic 

loads as presented herein are valid only in the normal flow region. 

The inculsion of the effects of the reversed flow region on the 
retreating blade can be incorporated in the model by first determining 
the region of reversed flow for the given flight condition. Re- 
calculating the airload expressions within this region using 

jz.fr ^ ~ C^3') 

^sr\c 

will model the influence of reverse flow. 

In Eqs. (39)-(42), the induced inflow field has not been 
specified. Typically it will have radial and azimuthal dependence. 
The high frequency and nonuniform nature of the induced Inflow 
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field can be incorporated using an analytical wake modeling program. 
On the other hand, if one is primarily concerned with low frequency 
variation, a dynamic Inflow model such as the one presented in Ref. 

20 may be used. 
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C. 3 Damping Loads 




E 


The blade structural damping loads are assumed to be 
of viscous nature (assumption 11). The distributed structural 
damping forces and moments written in the undeformed blade 
coordinate directions are 

-pl>3k,=^ 

respectively. 

C.4 Final Blade Equations of Motion 

The governing nonlinear partial differential equations 
of motion for the kth blade of an N-bladed rotor are obtained 
by substituting the inertia, aerodynamic, and damping loads 
into the equations of equilibrium, Eqs. (l)-(4). After some 
algebraic mainpulation and a final application of the ordering 
scheme, the equations for the lead-lag, flap and torsion 
deflections are obtained. These equations are presented in 
Appendix C. 

The boundary conditions which the blade deflections 
must satisfy depend on the type of rotor. For a hingeless 
rotor, the boundary conditions are 

(a) At the root (x^^^ = 0) ; 
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(b) At the tip (x^j^ " -O 

= /Wy = />^'g = T = Vy = Vg = O 

From Ref. 10, the blade cross section elastic moments and shears 
are given by; 

/riy ~ -BC - JTy y)«5-//7 9 q cos -f w, 

■" 'hXyyCGS 

+ Tx-fSin^c^ C4&) 

t B (X:^:^--Z‘yy') cS/n 9^0 05^^ 

- XXj c-^s Q-^ (4^) 

Vy ' - [e:(^ 

^ C T:2.:j.COS -f Xy y-3'/'1 ^ + 4>^>Xy^ 

- TyyCZosP^],^ “ -GXc^,^W,^ 
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f 


- S/f1 <^CoS 0(3 V, 4>,;^ + £ ( Xgj. ^//7 ’■0O 

■^XyyCOS <90 )w,^^4>,^ "■‘^3‘2 ^ 5 * 0 ) 


■== + 'JCs^ySIn 

+ 4>,^ cos(9'e, + G>y 4?,^ 

- [E(;i2X-TYY)'5'in0(S Ccps(9g Cv,^y V 4:>w^xx') 


+• E CXg£5*in‘^0(3 + lYV^^^^C5)Cw,xy-4^}yy'^ 

- Tx-f-^/n^Gl^y,'^ ^ ^>'<^3X "'^>X^3Z 

(S£> 


For an articulated rotor, 
blade tip remain the same, 
(a) At the root 


the natural boundary conditions at the 
The blade root boundary conditions are 
= 0): 


^r^v-w-O ; y»xy^>x) 


n^Y ~ 


where My and 


are given above. 


D. Fuselage Equations of Motion 

The governing fuselage equations of motion are determined by 
balancing rotor- and fuselage-induced forces and moments at the hub. 

These matching conditions between the N-bladed rotor and the fuselage ensure 
dynamic equilibrium of the overall system. 


Using the rigid body model of the fuselage and invoking D'Alembert's 
principle, the inertia force vector acting at the fuselage center of 
mass is given by 

/V ^ 

- - Mp + R-^ 0 ^ 2 ^ 

(Cr-:?:) 


Note that Eq, (52) is restricted (assumption 16) to flight conditions 
where the fuselage F-coordinate system is an inertial reference frame. 
The inertia moment vector, Q_, obtained from D'Alembert's principle and 
acting at the fuselage center of mass is given by 

~ 4l 

yu 

" - 1* ^ ~~ ^ ^ 

~^yQy 4 C^y'~J^x) 

4 9y ^2 ^ FX 

~ J Xy ^y 4- JZ~2X C 

-h 9^ - A< 

S-y C X-^ - Uy") *- 9x ^ j FV 
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-Xy/) + &x^Y ^ J^y"^ ^^y( ^X 

t$ 

i'Byi^y 3 ^y ^2 X.Py^jt')^ ^ FZ 

(r^3*> 


The fuselage-induced inertia force acting at the displaced hub is 
^H * fuselage-induced inertia moment acting at 

the hub, Q , is given by 

A/“ 

Q^.1 ^ Q f= f ^ 

xv ' /V ' Aj' 

~ “ f ^ ~ ^ ^2X ^zx C " ^ x^y 

~30y(9x') ^^y~ ^x'^ ^Y^z(-^y~^z^ 

i-ay(9^CTy-X^) + Mpr^^iRy ^R-^By) 

-f- M pr Xp C Ry ^ e. 

- [ Xy(9'y 4 -B-yBy ^ Bi-ex') 

+■ 6^^ <9x (■ 1-y ~ Xy) ■+■ By O-^ C X^2. " Xy') 

■+ BxBg.(y^~Xy) - M pypCR^ RyBy') 

~ /^ p ( Ry ~ RzBy)^ G. ^y 


31 



* e^e^(x^-ij.') ^ B-^ByCx^-iy') 

+ By&yCXx- Iy> ^MfX-fCRy -^y&jt') 

- Mptp ( Ry^-x + ^y'^Y^} ® P.X 

CS4') 

who."e, 1,'rom Eq. (10), 

>y- ’ cj'O 

The gravity loads acting at the hub, assuming the helicopter 
is in level horizontal flight, are given by 

^fiy 

- C Mp ^ Mp ) ^ COSOf^/2, ^ fz ^ 

4- A7^^<C c-o^^/z. Cxp - 
'^/'7 ^jz (."^F ^ ^pBy jzy 

Mp^s/n^^ C Xp9:z, + '^p&x') ^ 
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Gravity loads are considered to primarily influence the helicopter trim 
problem. The rotor blade gravity loads are considered to be negligible 
compared to the blade Inertia and aerodynamic loads and have been 
neglected in Sec. III. C. 

Other forces and moments acting on the fuselage and Induced at the hub 
may be expressed in general form ns 

(^T * ^X^RX ^Y'^RY c^’=t) 

These loads can represent airframe aerodynamic forces, including 

aerodynamie interference effects and tail rotor influence, or the effect 

of landing skids when in ground contact. This latter model could allow 

for the possibility of using this math model to investigate the dynamic 

21 

phenomenon of ground resonance. 

The final ingredient of the matching conditions are the rotor forces 
and moments induced at the hub, P and Q , respectively. The distributed 
loads acting on the kth blade expressed in the undeformed 3k-coordinate 
system are given by Eqs. (5) and (6). The individual components of these 
vectors were derived and presented in Secs. 111. C.l - III. C.3. 

Rewriting these loads in the rotating 2k hub plane system, p_, , and q , 
the three forces and moments Induced at the blade root of the kth blade 
are 
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C6:>0') 


^2K “ y Tz-k <=^'kk. 

0 ^ 

Szk ■ feu + ‘^^HkX Tzy] <^>^k 


where 


^Nk ““ ^ ^2Xk ■*■ ><w‘^3-xic •♦■ ^k*^3^k 


c^zy 

Transforming the forces and moments into the nonrotating rotor coordinate 

system, P^, » and Q , and performing a summation over all N blades of the 
^Rk ^Rk 

rotor yields the final expressions for the rotor-induced forces and moments 
at the hub. 







C^3y 

C ^4"') 


The final fuselage equations of motion are given by writing force 
and moment balance at the hub due to fuselage- and rotor-induced loads. 


Vf f P;z O 

Ay 

+■ Qq + Qf + (Xo =• O 

AJ " Ay ^ fj ' /V 


(.(oS-) 

C(p(o') 
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These equations, even after application of the ordering schfeSne, are very 
lengthy and reflect the complexity of modeling coupled rotor/fuselage 
dynamics. The equations themselves are given in Appendix D, 

IV. ROTOR /FUSELAGE SYSTEM 
EQUATIONS OF MOTION 

The final governing equations of motion include the three equations 
governing lead-lag, flap, and torsion deflections for each blade of the 
N-bladed rotor, Eqs. (C-2)-(C-4). Conseqvaently there are a total of 3N 
rotor equations of motion. Fuselage motion is governed by Eqs. (D-1) - (D-6) . 
Note that in the fuselage equations of motion, Eq. (11) was used to express 
the fuselage translational degrees of freedom in terms oi hub translations 
and rotations. 

f'x - &:£) 

-f- X jc <9^ + Syr* C 

^2 “ Tg - Xfr9-y 

Therefore the system motions are r^^, r^, r^, 0^, 6^, 6^, v^^, and 
k = 1, N. 

As written, the final equations are in second order partial differ- 
ential form. The spatial dependence may be removed by expressing the 
blade deflections as summations of elastic mode shapes and, by Galerkin's 


C&T) 

C(c>S>) 
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method, integrating the rotor equations to eliminate the spatial 
dependence. The selection of blade mode shapes will depend on the 
rotor model, hingeless or articulated. The number and type of mode 
shapes, used in each deflection representation, will directly effect 
the size of the final system of equations and the ability to determine 
system aeroelastic behavior. 

Several comments can be made about the nature of these second order 
nonlinear periodic coefficient equations. As derived, the rotor equa- 
tions of motion are written in a rotating reference system whereas the 
matching conditions between the rotor and the fuselage are written in 
a nonrotating system. Although this may cause some difficulty when 
visualizing results, this procedure allows for greater physical insight 
into the physics of the problem. 

It is possible to rewrite the rotor equations of motion in terms 

22 

of nonrotating rotor degrees of freedom using multiblade coordinates 
This transformation will replace some of the periodic coefficient terms 
in the final equations with constant coefficient terms. For some cases 
(hover, N ^ 3) this procedure will remove all periodic coefficient terms 
from the rotor equations of motion. However this is not the case for a 
helicopter in forward flight. Thus to be able to model an N-bladed rotor 
in forward flight, the rotor equations of motion will include periodic 
coefficients. Even after applying a multiblade coordinate transformation, 
it is possible to use a constant coefficient approximation for the remaining 
periodic coefficients. Although this approximation is good for many practi- 
cal cases, it is sensitive to advance ratio, operating state, and system 

23 

component frequencies . 
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The equations of motion presented herein can be used to study both 
aeroelastlc stability and response of coupled rotor and fuselage systems 
In hover or forward flight. There Is no restriction on the number of 
rotor blades. Additionally, due to the consistent formulation of the 
complete problem, these equations can be used to represent single isolated 
rotor blade behavior. By suppressing the hub motions in Eqs. (C-l)-(C-4), 
governing lead~lag, flap, and torsion equations of motion are obtained. 


V. CONCLUSIONS 

This paper presents the formulation of the governing equations of 
motion of a coupled rotor and fuselage helicopter math model. The 
derivation procedure yields a consistently derived set of second order 
nonlinear periodic coefficient partial differential equations. These 
equations may be used to study coupled rotor/fuselage dynamics in hover 
or forward flight. 
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TABLE 1 




DESCRIPTION OF COORDINATE SYSTEMS 


SYSTEM 

F 

R 

1 

2k 

3k 

4k 

5k 

V 


DESCRIPTION 

Fuselage axes. Origin fixed at undisturbed 
fuselage center of mass 

Rotor axes. Nonrotating. Origin fixed at 
undisturbed hub location. Directions 
coincident with F-system. 

Nonrotating rotor system with origin fixed 
in hub. Directions coincident with deformed 
R-system directions. 

Rotates with kth blade. Origin fixed in hub. 

Rotates with origin fixed in kth blade's 
pitch bearing. Preconed by • 

System yields undeformed blade coordinate 
directions. 

Origin fixed in deformed kth blade cross 
section at elastic center. System undergoes 
rigid body translation and rotation from 
3k-system. 

Origin fixed in deformed blade at elastic 
center. Rotated from 4k-system by removing 
elastic torsion rotation. 

Nonrotating. Origin fixed in hub. axis 

defines direction of oncoming air velocity 
due to forward flight of aircraft. 
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®2Zk 



3k * OAK 


^2k' ®2Xk 


Fig. 4 Rotor blade coordinate systems 



Fig. 5 Blade cross section configuration 
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o> 


^4k'®4Zk 



Fig. 6 Undef armed and deformed kth 
blade cross section 



Fig. 7 Deformed blade coordinate system 
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'^BZk 


Fig. 8 Airfoil velocities and steady 
aerodynamic forces on deformed 
kth blade cross section 


LuSk 



Fig. 9 Airfoil configuration for unsteady 
airloads calculation 
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APIPENDIX A . 


The Euler transformation matrices used in the derivation of the 
equations motion are given below. For rotations of 0(e) and higher, 

f 

the small angle assumption is used (cos sin 0« 0 ). 



fe^vJ rcosHf. 0 few 


» ' 


COS^p 

0 

s'in/Sp 


• • 
^2Xk 

^ ^3Yk 

»• — 
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‘^ZYk «' 
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» < 
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^2Ek 


• • 

®4Vk 
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'^k’y 
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A 
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r 

^5Xk 


1 

0 
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^4Xk. 

“* ^5Yk 

r — 

0 
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^4Vk'' 

®5Zk 


0 
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^4Zk 
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- 


49 




(A. 4-) 


(A.5) 


C A. 6^ 




APPENDIX B 


From basic definitions, Fig. 5, and the assumptions of section 
II. B the following identities are used in the derivation. 
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*1 
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CB.|5) 

where 
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(B.\(i') 

4k = 



CB.17) 
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CB./8) 
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APPENDIX C 


The governing nonlinear partial differential equations of motion 
for kth blade of an N-bladed rotor are presented below. They represent 
equations of equilibrium in the axial direction and the lead-lag, flap, 
and torsional deflections, respectively. 

C :L t O Cc-l') 

^2:Z: CIOS'*- 4. XyySin'*^(9^^.':> 

+• E C - Xy y") 2 

0 

^ ^ ~ ■^>77 32^c='Os- ■> X 

C-OsVl^^^ <f-ry 4- 'X<Ny_Oy + 9^ 

+ /Sp Ce + y,^')') 4- 5-|n^|^(^ Ty + 2W,^<9y 

• • 

•<' c W|^ + /3p Ce + 2/SpWv, 

7^A3Vk “ ^ Cc-z:) 
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E (f. Ig^ sin' ^Gk •" ^VY '^Gk‘> '^k .X V 

■*■ ^ ( Xg g “ Xyy^ 4*V< '^U • w ^ ^Gk 

- C^ x-r 4- 

C >x ”^k”^ > X -O- j^ ^x ^ ^ t" S7/7 

^^G>k ^r^T3S -"x 

+ 777_CT2- Ccds *s ~ X V — Vj^ - 2. ^y 

+ y3p^y " C ^ ^y - ^y') C C i- 

+■ Sin - Z 0y^k ~ ^k'^Y -k /5pPY 

- (20y + (e+X^;^> -fg^ -Wi^ ~'2jSyC</u. 

- (^e-VXj^'>') -/Sp (e.-VYvJ> ~ ^x 

- ^SF^k-^ -pAS-^k " ^ 
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ir -TV? C.O-S C e,4 X\^ C Qy ~ 2 

- CevKO CZBy -V - r^ - <i7u 

- 2\>j^ C/Sf-^wv^,^^ 4 <i>^C Vi^-Vj^') Ce+xO C-/3p 

Vvo^<^w:^] [ V|^ 

4 Vi^,^ Ce 4 Xi^i] - -oJ^ Xp 4 
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The elastic axial degree of freedom is defined by 



Cc-s') 
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The aerodynamic quantities P.„„. » 

A JX K 

Eqs. (41) and (42), respectively. 
Eqs. (C-2) - (C-4) is given by 


"ASZk • >Xk 

The coefficient appearing in 



J' 7T7 C <2. 4- dx 


-SO- 


Cc-6) 


Appendix D 


The follovdng equations represent force and moment equilibrium at the 
rotor hub. 


"/W^f py - EpCl^y 4- 

x-pdg - &yO^- 2 dy&z-dY&^) 

Fy ^ "Pry ~ o (D. 2 ) 


~Mp C Tzt - Xr^y^ “ 

F-^ + pRr^ =■ 0 C^^-s") 


-(Z,-2T^^dy)§, ' r^xC^^-^y^y-3dy/yy) 
•f ('Xy -:T y)(^z(3v CX^-X^)&yd^ 

^ CZ^- Xy)&yd^ - -V- Xp(<^2-§-y0y'^ 

C<^x ■ ~ ^2 ^x] 

- A^y=: X^[ + py 0-y + + Oy^g) + 
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4 4- Qfz,y — O 

-Xy 0y - X^x ^^2 ^z “ ^x<^x - ^>T) ^ (^y-^y) 

■« 

■f ^Xy "Xg.") - •^z') ■*‘■'^'^^^[^^2 

- <^>^p- -f -^pO^y^dy -i-ry(9y] ^Mp-^plc^ 4- X p9s( 9y 

- Fg (^ <^Y ^ ^ ^2 ■'■ 

^ C ^ p’ ■” ^F" 4~ <s' 1 ^r* pj_ ^2^4- Vp- ^^y') j[ 

+ Qy ^?RV ' — C^.-^ ) 

-Xj&g -X2 x( 4 +ez®r +3(9yeg +?i9y4Vl9yi9/-^2’^x) 

+ Cry-:z2^&x©y <- (Xy-r^)4'^Y ~ 

+ x^Oa V 4 -2i9yi9j - ey^2)]4-Mp*p[rx6x 
+ ry&Y + XpSaOv* ZpC-SyBx ’■®Y®a^] 

+ f^F^s\oo('e, ( Xf^ 9" 2 ^ 4 + ©2 ■'■ Qr.2. ~ ^ 
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4 


The generalized fuselage forces and moments, F^, Q^, j • X, Y, Z, 

are left undefined# The rotor-induced forces and moments, P„ Q.„ , 

Rj, Rj 

j » X, Y, Z, are defined in Appendix E. 
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APPENDIX E 


The following are the rotor~induced forces and moments acting 
at the hub in the R-systero directions. 

N 

o *■ 

• • 

*♦■ 2P-^) + ^»~Y C3/^p-^^k.x‘> 

• • » * - _■ . t 1 

4 ■'" '^■g- ■' ^ 

4- Vs/j^"^ - 2. ^/p<. Ce4-*^,^'^/3p^g J ■*■ cc>s4’^^s'in4^^j^^ (9yV|^,;^ 
4 :2^ ^2 Ak + ^ '^■fe ^ 

+ fy ( Ai,4 4 3/<y6pC v7j^ 4 Ce.4-xO^-^"^ 

4 ^ ( ^Ic ^ ^ e-V-kj^l ( 1. + (9-y 

- C e- *v- X yj) V Y +• 2r~Y|c3pCe-4-k J 

4 CO s 4’|^ [^ - ^ pCe4X)<r)(^2V|^ + C'e-4x^)C’i42 ^ •+■ 
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<• 3in+u[ -3^0y c% ¥ ca.^.so'> +'j^, CAu+n +^1)1,'^ 

” '^k^x'^k>x -V r^ + ^j^') 

+ Ce -4-xO ^ 

+• Ce4>x.^'^©Y C + CS4-Xi^‘)^ 

o- < S;ri [ -^ ,: b ^ 5J^-^J 

<■ S/n’-^t ‘^“'■‘^'k [, -/^/Sp&Y C b -S^'^'] 

+ cos’’4'k[ 2/SpCy^W;^,^ - CSvx^&y 
+ i:^i-J^k'>^x -><-CG>-c^k+4>0') -/Speyfe+x^^-^Ci-Sj^-)] 

. ->^/Sp ^"ic>. + /Sp'^Ci 

->‘-<^^Sk+'i|x'^C|;W-Xp,')^p + CeH-X^,-^(9Y^pi 

- C e + X „> e-y ^ p (|. . x^-)] p Cw„^-v^p-)(r£ 

■*■ XjJ) C B - x^') ~'(^f C^Jt~‘^^(^^ b“X^ 

- ^ /Sp C^^i&y - A^. - Ts. - +^ e-y I 4,^. 

a b-xO C >^k,« ^/3p 1 6 ^;,+ 4>^A] 

+ (Py f e4-X|^')fk/^,^ +/SpK£-Xp') + 6 *y fet-xt.A'^esk'ni'tACl^ 

“ ><a') ~ ^ Xk + ^kf^ a ^ '*' 
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+ XEOoj©g^ + 4P^©-J + +^u«9z 
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N 

^/ez ' 2 
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